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Abstract. We investigate the centrifugally driven curvature drift instability to study how field lines twist close to the light 
cylinder surface of an AGN, through which the free motion of AGN winds can be monitored. By studying the dynamics of 
the relativistic MHD flow close to the light cylinder surface, we derive and solve analytically the dispersion relation of the 
instability by applying a single particle approach based on the centrifugal acceleration. Considering the typical values of AGN 
winds, it is shown that the timescale of the curvature drift instability is far less than the accretion process timescale, indicating 
that the present instability is very efficient and might strongly influence processes in AGN plasmas. 
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INTRODUCTION 

For studying AGN winds the fundamental problem re- 
lates to the understanding of a question: how the plasma 
goes through the Light Cylinder Surface (LCS), which 
is the hypothetical zone where the linear velocity of ro- 
tation equals the speed of light. This implies that the 
plasma particles, which move along quasi-straight mag- 
netic field lines in the nearby area of the LCS, must reach 
the speed of light. Generally speaking no physical sys- 
tem can maintain such a motion and a certain twisting 
process of the magnetic field lines must operate on the 
LCS. On the other hand if the trajectories are given by 
the Archimedes spiral, then the particles can cross the 
LCS avoiding the light cylinder problem yj]. An addi- 
tional step in this investigation is to identify the appropri- 
ate mechanism that provides the twisting of the magnetic 
field lines, giving rise to the shape of the Archimedes 
spiral, and in turn insures that the dynamics is force-free. 

Since the innermost region of AGNs rotates, the role of 
the Centrifugal Force (CF) appears interesting. The cen- 
trifugally driven outflows have been extensively studied. 
Generalizing the work developed in [2] it was shown that 
due to the centrifugal acceleration, electrons gain very 
high energies with Lorentz factors up to y ~ 10 8 JH0]. 
This implies that the energy budget in the AGN winds is 
very high. 

The centrifugally driven parametric instability was 
first introduced in yllg] for the Crab pulsar and AGN jets 
respectively. Another kind of the instability which might 
be induced by the CF is the so called Curvature Drift 
Instability (CDI). In |7[] the two-component relativistic 
plasma has been considered to study the role of the cen- 
trifugal acceleration in the curvature drift instability for 
pulsar magnetospheres. The investigation demonstrated 
high efficiency of the CDI. To investigate the twisting 




FIGURE 1. Two orthonormal bases are considered: i) cylin- 
drical components of unit vectors, (e$,efi, e v ); ii) unit vectors 
of the system rigidly fixed on each point of the curve, (e r ,eg, 
e v ), respectively. C is the center of the curvature. Hereafter, this 
set of coordinates is referred to as the field line coordinates. 



process of magnetic field lines due to the CDI, we apply 
the method developed in 10, Ht 1 to AGN winds. 

The paper is arranged as follows. In Sect. 2, we intro- 
duce the curvature drift waves and derive the dispersion 
relation. In Sect. 3, the results for typical AGNs are pre- 
sented and, in Sect. 4, we summarize our results. 



MAIN CONSIDERATION 

We begin our investigation by considering the two- 
component plasma consisting of the relativistic electrons 
with the Lorentz factor y e ~ 10 5 ~ 8 (see d 01) and the 
bulk component (protons) with Yb ~ 10. Since we are 
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FIGURE 2. The dependence of logarithm of the instability 
timescale on the normalized wave length. The set of parame- 
ters is YeO = {10 5 ; 10 6 ; 10 7 ; 10 8 }, R B w Ri c , n e0 = 0.001cm -3 , 
= l00R lc , L/L E = 0.01, L E = I0 46 erg/s is the Eddington 
luminosity for the given AGN mass. 



FIGURE 3. The dependence of logarithm of the instability 
timescale on the density normalized by the medium density. 
The set of parameters is y e0 = {10 5 ; 10 6 ; 10 7 ; 10 8 }, R B ~ Ric, 
X e Q — R[ c , Xa = lOOR/c and L/Le = 0.01. Here, n e Q is normal- 
ized by the intergalactic medium density, n m ~ 1cm -3 



interested in the twisting process, we suppose that ini- 
tially the field lines are almost rectilinear to study how 
this configuration changes in time. 

We express the equation of motion in the cylindrical 
coordinates (see Fig. [TJ and we start by considering the 
Euler equation: 

% + (v a V)p a = -c 2 7a ^ + SB. (e + iv„ x bV 
at m a \ c J 

(1) 

the continuity equation: 

^+V(n a v a )=0, (2) 
and the induction equation: 

VxB=I§ + ^£j a , (3) 

C C a=e,b 

where a — {e,b}, | = \J\ — Q. 2 R 2 /c 2 „ p a is the mo- 
mentum, v a - the velocity and y a is the Lorentz factor of 
the relativistic particles. 

We linearize the system of equations Eqs. (Q][3]), per- 
turbing all physical quantities around the leading state 
xj/ w ip0 + ^i ; where xj, = {„ jVjPj e,B} and ¥ l {t,r) °c 
^(Oexp^kr)]. 

Then, from Eqs. (JTJ3J we derive the linearized set of 
equations governing the CDI: 

^-i{k x u a + k^v y )/4=— * X, (4) 
at 11 m a 11 



- i(k x ii a + k^v^ )n l a = ik x n° a v l ax , (5) 
-ik$cB l r =4n £ qa(naV l ax + n l a u a ). (6) 

a=e,b 

By u a — 7(j^v?/((Db2?b)> we denote the curvature drift 
velocity along the x axis, where co a B = QaBo jm a c, Rb 
is the curvature radius of magnetic field lines, v., « 

ccos(iif), Bq = \j2Lj (Ri c c 2 ), L is the luminosity of the 
AGN and R[ c = c/Q. - the light cylinder radius. In de- 
riving Eqs. dU|6]l, the wave propagating almost perpen- 
dicular to the equatorial plane, was considered and the 
expression v\ « cE\ /Bq was taken into account. For sim- 
plicity, the set of equations are given in terms of the co- 
ordinates of the field line (see Fig. |T}. 

We express v l ax and n l a in the following way: 

V^Va****®, (7) 

n l a =N a e^\ (8) 

where 

A a ,(0 = ^ + ^sm(2^), (9) 

A^(f) = ^sin(nf). (10) 

Then, by substituting Eqs. (jTJ and ((HJ into Eqs. (@]|6]i, 
it becomes straightforward to solve the system for the 
toroidal component and find a corresponding increment 
of the instability: 
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FIGURE 4. The dependence of logarithm of the instabil- 
ity time scale on log(y e o)- The set of parameters is Rr ~ 
Rio "co/«m = 0.001, XeQ = R[ c , = 100i?/ c and L/L E = 
{0.01;0.05;0.3;1}. 
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FIGURE 5. The dependence of logarithm of the instability 
time scale on log(y e o). The set of parameters is Rb ~ Ri c , 
neo/n m = {0.001;0.01;0.1;1}, A e0 = Ri c , = 100J?/ C and 
L/L E =0.01. 



RESULTS 

We investigate the CDI growth rate in terms of the wave- 
length, the density of relativistic electrons, their Lorentz 
factors and the AGN bolometric luminosity. 

In studying the behaviour of the instability as a func- 
tion of the wavelength, we examine the typical AGN pa- 
rameters: M B h = 10 8 x M , £2 = 5 x lCT 5 ^ 1 and L = 
10 44 erg/s, where Mbh is the AGN mass, M Q is the solar 
mass and L is the bolometric luminosity of the AGN. 

We consider Eq. (fTTT i and plot the logarithm of the in- 
stability timescale t = 1/T versus the wavelength. The 
present consideration is based on the centrifugal accel- 
eration. As shown in [3], due to the CF, the relativistic 
particles can reach very high Lorentz factors. For this 
purpose, it is reasonable to investigate the efficiency of 
the instability in terms of the wavelength but for differ- 
ent values of Lorentz factors. Fig.|2]shows the mentioned 
behaviour for different parameters. Different curves cor- 
respond to different values of Lorentz factors. For the 
given range of X x and different values of 7,0, the CDI 
time scale varies from ~ 10 7 s (X x /Ri c = 1, Je0 = 10 8 ) to 
~ 10 9 (X x /R lc = 2, 7e0 = 10 5 ). 

In Fig. [3] the plots of log(x) versus the AGN wind 
density illustrate that the timescale is a continuously 
decreasing function of n e o jn m . As we see from the figure, 



X varies from ~ 10 9 s (n e o/n m = 0.001, y e o = 10 5 ) to 

- 10 6 s (n e0 /n m = 1, 7,0 = 10 8 ). 

In Fig.|U the plot of log(x) versus log(y e o) is shown 
for different luminosities. As we see, the instability 
timescale varies from ~ 10 10 .? (y e o = 10 5 , L/Le = 1) to 

- 10 8 .? (7 e0 = 10 7 , L/L E = 0.01). On the other hand, the 
plots for different luminosities illustrate another property 
of T: by increasing the luminosity of the AGN, the corre- 
sponding instability becomes less efficient. 

To observe this particular feature more clearly, we 
consider how in Fig. the dependence of log(x) on 
L/Le is clearly evident for different values of densities. 
From the plots, it is seen that by increasing the luminos- 
ity, the timescale continuously increases. For the afore 
mentioned area of quantities, the timescale varies from 

- 10 7 s (L/L E =0.01, n e0 /nj= 1) to - 10 10 s (L/L E = 1, 
n e o/nj = 0.001). 

We observe from the present investigation that the 
instability timescale varies in the following range: x G 
{10 6 ;10 10 }s. To specify how efficient the CDI is, it 
is pertinent to examine an accretion process, estimate 
its corresponding evolution timescale, and compare this 
value with that of the CDI. 

Considering the problem of fuelling AGNs 
H it was showed that the accretion timescale 
can be estimated by the following form t evo i = 

3.53 x lO l3 {lOL/L E y 22/21 (lO-*M BH /M Q y 4/21 s. 
As is clear from this formula, the accretion evolution 
timescale depends on two major AGN parameters, 
the luminosity (L) and the AGN mass (Mbh)- There- 
fore, it is reasonable to investigate t ew ,i versus L 
and Mbh- Let us examine the following ranges of 



variables: M 9 = M BH /(10 9 x M ) e {0.01; 1} and 
L/Lg S {0.01; 1}. Then it is easy to show that the 
minimum value of the evolution timescale is of the 
order of ~ 10 l2 s, which corresponds to M9 = 1 and 
L/Le = 1, whereas the maximum value, approximately 
~ 10 15 s, corresponds to the following pair of variables 
M 9 =0.001 and L/L E =0.01. 

As has been found, % varies in the range ~ (10 6 — 
10 10 )s, whereas the sensible area of t evo i is ~ (10 12 — 
10 15 )s. Therefore, the instability timescale is less than 
the evolution timescale of the accretion by many orders 
of magnitude, which implies that the linear stage of the 
CDI is extremely efficient. 

The twisting process of magnetic field lines requires 
a certain amount of energy and it is natural to study 
the energy budget of this process. For this reason we 
have to introduce the maximum of the possible luminos- 
ity L max = Mc 2 and compare this with the "luminosity" 
corresponding to the reconstruction of the magnetic field 
configuration L m = AE m /At ~ AE m / x, where AE m is the 
variation in the magnetic energy due to the curvature drift 
instability. 

We consider a AGN of the luminosity, L = 10 45 erg /s, 
then, the accretion provides the following maximum 
value: L max — L/0.1 = 10 46 erg/s. On the other hand, if 
the process of sweepback is realistic, the magnetic "lu- 
minosity" cannot exceed L max . The magnetic "luminos- 
ity" can be expressed by following L m = B 2 .R? 1c k/{Ati%), 
with B, = B^exp(t/l), where B® is the initial pertur- 
bation of the toroidal component of magnetic field and 
K = AR/Ri c << 1 represents the non-dimensional thick- 
ness of a thin spatial layer close to the LCS. 

We introduce the initial non-dimensional perturbation, 
8, defined to be 8 = B^/Bq, where by Bq we denote 
the induction of the magnetic field in the leading state. 
By considering the following set of parameters y e Q = 
{10 5 ;10 6 ;10 7 ;10 8 }, R B « R lc , n M = 0.001cm~ 3 , X$ = 
lOORic, X x = Ric and L = 10 45 erg/s, we investigate the 
behaviour of L m /L max versus the initial perturbation for 
the characteristic timescale it w t). One can see that, 
L m /L max varies from - (8 = 0) to - 2.3 x 10~ 9 (8 = 
0.1, Yeo = 10 s ). Therefore, L max >> L m which means 
that only a tiny fraction of the total energy goes to the 
sweepback, making this process feasible. 

SUMMARY 

We summarize the principal steps and conclusions of our 
study to be: 

1 . Considering the relativistic two-component plasma 
for AGN winds, the centrifugally driven curvature 
drift instability has been studied. 



2. Taking into account a quasi single approach for the 
particle dynamics, we linearized the Euler conti- 
nuity and induction equations. The dispersion rela- 
tion characterizing the parametric instability of the 
toroidal component of the magnetic field has been 
derived. 

3. By considering the proper frequency of the curva- 
ture drift modes, the corresponding expression of 
the instability increment has been obtained for the 
light cylinder region. 

4. The efficiency of the CDI has been investigated 
by adopting four physical parameters, namely: the 
wavelength, flow density and Lorentz factors of 
electrons, and the luminosity of AGNs. 

5. By considering the evolution process of accretion, 
the corresponding timescale has been estimated for 
a physically reasonable area in the parametric space 
L — Mbh- It was shown that the instability timescale 
was lower by many orders of magnitude than the 
evolution timescale, indicating extremely high effi- 
ciency of the CDI. 

6. Examining the instability from the point of view of 
the energy budget, we have seen that the sweepback 
of the magnetic field lines requires only a small 
fraction of the total energy, which means that the 
CDI is a realistic process. 
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